Effects of static charging and exfoliation of layered crystals 
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Using first-principle plane wave method we investigate the efTects of static charging on structural, 
electronic and magnetic properties of suspended, single layer graphene, graphane, fluorographene, 
BN and M0S2 in honeycomb structure. The limitations of periodic boundary conditions in the 
treatment of negatively charged layers are clarified. Upon positive charging the band gaps between 
the conduction and valence bands increase, but the single layer nanostructures become metallic 
owing to the Fermi level dipping below the maximum of valence band. Moreover, their bond lengths 
increase leading to phonon softening. As a result, the frequencies of Raman active modes are 
lowered. High level of positive charging leads to structural instabilities in single layer nanostructures, 
since their specific phonon modes attain imaginary frequencies. Similarly, excess positive charge is 
accumulated at the outermost layers of metallized BN and M0S2 sheets comprising a few layers. 
Once the charging exceeds a threshold value the outermost layers are exfoliated. Charge relocation 
and repulsive force generation are in compliance with classical theories. 



I. INTRODUCTION 

Single layer graphene, [T] graphane CH,[5J [3] fluoro- 
graphene CF,[1[S] BN[il3 and MoSzlHlin] have displayed 
unusual chemical and physical properties for future nan- 
otechnology applications. Furthermore, the properties of 
these nanomaterials can be modified by creating excess 
electrical charge. For example, linear crossing of bands 
of graphene at the Fermi level gives rise to electron-hole 
symmetry, whereby under bias voltage the charge carriers 
can be tuned continuously between electrons and holes 
in significant concentrations. [TU] This way, the conduc- 
tivity of graphene can be monitored. Similar situation 
leading to excess electrons or holes can also be achieved 
through doping with foreign atoms. pTMl5) Layered ma- 
terials can be exfoliated under excessive charging, which 
is created by photoexcitation for very short time. |T6l [17] 
It is proposed that the femtosecond laser pulses rapidly 
generate hot electron gas, which spills out leaving behind 
a positively charged graphite slab. Eventually, charged 
outermost layers of graphite are exfoliated. [17] 

Recently, the effects of charging of graphene have been 
treated in different studies. Ekiz et [THj showed that 
oxidized graphene domains, which become insulator upon 
oxidation, change back to the metallic state using elec- 
trical stimulation. Theoretically, based on the first prin- 
ciples calculations, it has been shown that the binding 
energy and magnetic moments of adatoms adsorbed to 
graphene can be modified through static charging. [151 [T9] 
Possibility of transforming the electronic structure of one 
species to another through gating modeled by charg- 
ing has been pointed out. [20] It is argued that diffusion 
of adsorbed oxygen atoms on graphene can be modi- 
fied through charging. [3T] We found that pseudopotential 
plane wave calculations of charged surfaces using periodi- 
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cally repeating layers are sensitive to the vacuum spacing 
between adjacent cells and have limited applicability. [15] 
In this paper we investigate the effect of static charging 
on suspended (or free standing) single layer nanostruc- 
tures, such as graphene, graphane (CH), fluorographene 
(fully fluorinated graphene) (CF), boron nitride (BN) 
and molybdenum disulfide (M0S2). All these honeycomb 
nanostructures have two dimensional (2D) hexagonal lat- 
tice. First, we examine how the size of the "vacuum" po- 
tential between layers affects the calculated properties of 
the negatively charged single-layer nanostructures when 
treated using periodic boundary conditions. We then in- 
vestigate the effect of charging on the electronic energy 
band structure and atomic structure. We show that the 
bond lengths and hence 2D lattice constants increase as 
a result of electron removal from the single layer. Con- 
sequently, phonons soften and the frequencies of Raman 
active modes are lowered. As a result of electron removal, 
three-layer, wide band gap BN and M0S2 sheets are met- 
allized and excess positive charge is accumulated mainly 
at the outermost atomic layers. Owing to Coulomb force 
those layers start to repel each other. When exceeded 
the weak van der Waals (vdW) attraction, the repulsive 
Coulomb force initiates the exfoliation. 



II. METHOD 

The present results are obtained by performing 
first-principles plane wave calculations carried out 
within spin-polarized and spin-unpolarized density func- 
tional theory (DFT) using projector-augmented wave 
potentials. 22 The exchange correlation potential is ap- 
proximated by Generalized Gradient Approximation. [23] 
For a better account of weak interlayer attraction in lay- 
ered crystals, van der Waals (vdW) interaction is also 
taken into account. [24) A plane- wave basis set with ki- 
netic energy cutoff of 500 eV is used. All atomic posi- 
tions and lattice constants are optimized by using the 
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conjugate gradient method, where the total energy and 
atomic forces are minimized. The convergence for energy 
is chosen as 10^^ eV between two steps, and the maxi- 
mum force allowed on each atom is less than 0.01 eV/A. 
The Brillouin zone (BZ) is sampled by (15x15x5) special 
k-points for primitive unit cell. Calculations for neutral, 
as well as charged systems are carried out by using VASP 
package.!^ 

Two-dimensional single layers or slabs and a vacuum 
space s between them are repeated periodically along the 
perpendicular z-direction. The amount of charging is 
specified as either positive charging, i.e. electron deple- 
tion {Q > 0), or negative charging, i.e. excess electrons 
{Q < 0), in units of ± electron (e) per unit cell. Av- 
erage surface charge density is specified as a = Q/A, 
i.e the charge per unit area, A, being the area of the 
unitcell. Normally, periodic boundary conditions real- 
ized by repeating charged supercells has a divergent elec- 
tric potential energy and has drawbacks and limitations, 
which have been the subject matter of several studies 
in the past. To achieve the convergence of electronic 
potential, additional neutralizing background charge is 
applied. [221 [37] Recently, error bars in computations due 
to compensating charge have been estimated. [20j The 
dipole corrections can be carried out for cubic struc- 
tures, if a finite electric dipole moment builds in the 
unit cell.[5Hllin] Monopole and dipole corrections are also 
treated self-consistently. [30] Various charged structures 
have been also treated by using different approaches and 
computational methods. |31f[57) Owing to those theoreti- 
cal advances, studies on charged systems can now reveal 
useful information, when treated carefully. 



III. CHARGING OF SUSPENDED SINGLE 
LAYERS 

The negative and positive charging of suspended sin- 
gle layer graphene, CH, CF, BN and M0S2 are treated 
using supercell method. In Fig. [l] (a) we describe M0S2 
single layers, which are periodically repeated along the z- 
direction and separated by a vacuum spacing s between 
the adjacent outermost sulfur planes. In Fig. [T] (b) and 
(c) the self-consistent electronic potential energy, Vei{r) 
is averaged in the planes perpendicular to the z-axis to 
obtain planarly averaged ID potential energy Vei{z) for 
different values of s. 

In the vacuum region, the electronic potential energy 
Veiiz) strongly depends on the vacuum spacing s. For an 
infinitely large single plane having excess charge Q > 
per cell, the potential energy in the vacuum region is 
linear, if it is not periodically repeating. Thus, as z ^ 00, 
Vei{z — 00) — >■ -l-oo as schematically shown in Fig. [l] 
(b). However, for a periodically repeating single layers 
(within the periodic boundary conditions) the potential 
energy is symmetric with respect to the center of vacuum 
spacing and it passes through a maximum at s/2. The 
maximum value of the potential increases with increasing 




z(A) 

Figure 1; (Color online) (a) Description of supercell geome- 
try used to treat 2D single layer M0S2. c and s are supercell 
lattice constant and vacuum spacing between adjacent layers. 
The z-axis is perpendicular to the layers, (b) Self-consistent 
potential energy of positively charged (Q > per cell), pe- 
riodically repeating M0S2 single layers, which is planarly av- 
eraged along z-direction. Vei{z) is calculated using different 
vacuum spacings s as specified by inset. The planarly aver- 
aged potential energy of a single and infinite M0S2 layer is 
schematically shown by linear dashed lines in the vacuum re- 
gion. The zero of energy is set at the Fermi level indicated 
by dash-dotted lines, (c) Vsi{z) of negatively charged (Q < 
per cell) and periodically repeating M0S2 single layers. Av- 
eraged potential energy of infinite M0S2 single layer is shown 
by linear and dashed line in the vacuum region, (d) Variation 
of Vei{z) and total potential energy including electronic and 
exchange-correlation potential, Vei{z) + Vxc{z), between two 
negatively charged M0S2 layers corresponding to Q=-0.110 
e/unitcell before the spilling of electrons into vacuum. The 
spacing between M0S2 layers is s = 20 A. (e) Same as (d) but 
Q=-0.114 e/unitcell, where the total potential energy dips 
below Ep and hence excess electrons start to fill the states 
localized in the potential well between two M0S2 layers, (f) 
Corresponding planarly averaged charge density A. Accumu- 
lation of the charge at the center of s is resolved in a fine 
scale. Arrows indicate the extremum points of Vei{z) in the 
vacuum region for Q > and Q < cases. 
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sin Fig.igb). 

In contrast, for a negatively {Q < per cell) charged 
and infinite M0S2 single layer, a reverse situation occurs 
as shown in Fig. [T](c). Namely Vei{z — 00) — —00 lin- 
early, if M0S2 single layer is not periodically repeating. 
Notably, the energy of a finite size, single layer nanos- 
tructure (i.e. a flake) does not diverge, but has finite 
value for large z both for Q > and Q < cases. On 
the other hand, for periodically repeating single layers 
within the periodic boundary conditions, potential ener- 
gies are symmetric with respect to the center of vacuum 
spacing and they passes through a minimum at s/2. This 
way a potential well is formed in the vacuum region be- 
tween two adjacent layers. Normally, the depth of this 
well increases with increasing negative charging and s. 
At a critical value of negative charge, the self-consistent 
potential energy y(r) including electronic and exchange- 
correlation potential energies dip below the Fermi level 
(even if Vei{z) > Ep) and eventually electrons start to 
occupy the states localized in the quantum well. Such a 
situation is described in Fig. [T] (d)-(f). Of course, this 
situation is an artifact of the method using plane wave 
basis set and the repeating layers separated by the vac- 
uum space s. Despite that, the method may provide 
reasonable description of the negatively charged layers 
until the minimum of the well dips below the Fermi level. 
According to this picture, the escaping of electrons out 
of the material is delayed in relatively short s. On the 
other hand, the interaction between layers prevents one 
from using too short s. Earlier, this limitation of the 
method is usually overlooked. The critical value of neg- 
ative charge depends on s value. It should be noted that 
for s=20 A, electrons start to escape from the graphene 
layer for (5=-4. 03x10^'^ e/cm^, even though larger doping 
of 4x10^* e/cm^ has been attained for garphene on Si02 
substrate. [55] 

In the case of positive charging, even if Vei{z) is not lin- 
ear and does not increase to -l-oo, the periodic boundary 
conditions using sufhciently large s can provide a realis- 
tic description of charged systems, since the wave func- 
tions in the vacuum region rapidly decay under high and 
wide potential barrier. Therefore, the calculated wave 
functions and electronic energies are not affected even 
if Vei{z) is smaller than the electronic potential corre- 
sponding to infinite vacuum spacing. We demonstrate 
our point of view by solving directly the Schrodinger 
equation to obtain the wave functions and energy eigen- 
values for the planarly averaged ID potentials of sin- 
gle layer and 3-laver graphene corresponding to s=12.5, 
25, 50 A in Fig. [2| One sees that the large difference, 
AVei(z) ■'(z)-V '(z) do not affect the 

occupied states at their tail region in the vacuum spac- 
ing; the energy difference is only 5 meV (which cannot 
be resolved from the figure) between smallest and largest 
vacuum spacing s, which is smaller than the accuracy 
limit of DFT calculations. As one expects, the depen- 
dence on the vacuum spacing increases for excited states, 
which have relatively larger extension and hence they are 
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Figure 2; (Color Online) (a) Energy eigenvalues of the oc- 
cupied electronic states, Ei and corresponding [^^^(z:)^ are 
obtained by the numerical solution of the Schrodinger equa- 
tion of a planarly averaged, ID electronic potential energy of 
single layer graphene for s=12.5 A, 25 A and 50 A shown by 
dashed lines, (b) Same as (a) for 3-layer graphene. Zeros of 
[^^^(z)!^ at large z are aligned with the corresponding energy 
eigenvalues. 

affected from AVei{z). 

By taking the above limitations of the method in neg- 
ative charging into account, we now examine the effect 
of charging of single layers of graphene, CF, CH, BN and 
M0S2 on their electronic structure and bond lengths. In 
Fig. [3] the changes in band structure with charging are 
significant within DFT. For example, the band gap (i.e. 
the energy gap between the top of the valence band and 
the minimum of the conduction band) of neutral single 
layer BN increases from 4.61 eV to 5.12 and to 5.54 eV as 
Q increases from Q=0 to -1-0.2 e/cell and to -t-0.4 e/cell, 
respectively. The increase of the band gap occurs due 
to the fact that the electronic potential energy becomes 
deeper with increasing electron depletion. For Q > 0, the 
Fermi level dips in the valance band and creates holes. 

In contrast, parabolic free electron like bands, which 
occur above the vacuum level in the neutral case, start 
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Figure 3: (Color Online) Energy band structures of 2D single 
layer of graphene C, fiuorographene CF, graphane CH, BN 
and M0S2 calculated for Q — +0.2 e/cell, Q = (neutral) 
and Q = —0.10 e/cell. Zero of energy is set at the Fermi level 
indicated by dash-dotted lines. The band gap is shaded. Note 
that band gap increases under positive charging. Parabolic 
bands descending and touching the Fermi level for Q < are 
free electron like bands. Band calculations are carried out for 
s=20 A. 



a wide tunneling barrier, even if Vei{z) is lowered below 
the Fermi level in vacuum for large z. 

Incidentally, for both Q > and Q < 0, the spin- 
polarized calculations carried out for single layers of 
graphene, CH, CF, BN and M0S2 did not yield any mag- 
netic state as a result of charging. 

Another crucial effect of charging appears in the vari- 
ation of the bond lengths with Q. As shown in Fig. |4] 
(a) the bond length or lattice constants of single layer 
graphene, BN, CH, CF and M0S2 increase with increas- 
ing positive charge density a. The elongation of the bond 
length is slow for small a, but increases quickly when 
CT > 1 C/m^. The bonds get longer when the electrons 
are removed from the system and hence bonds become 
weaker. The contour plots of total charge density in a 
plane passing through C-C and B-N bonds of graphene 
and BN honeycomb structures in Fig.|4](b) and (c), show 
that the charge density between atoms becomes weaker 
with increasing electron depletion. Weakening of bonds 
can have crucial consequences as phonon softening and 
is observable through Raman spectrum. In fact, the Ra- 
man active mode of graphene calculated by using density 
functional perturbation theory is at 1524 cm~^ and shifts 
down to 1510 cm~^ for Q=+Q.2 e/cell, and to 1311 cm~^ 
for (5=0.4 e/cell. To confirm whether the elongation of 
bonds dominates the Raman shift of graphene, we calcu- 
lated the Raman active modes of neutral graphene hav- 
ing the same lattice constant of graphene when charged 
by Q~~\-QA e/cell. In this case the Raman active mode 
shifted to 1274 cm~^, which is close to the Raman ac- 
tive mode calculated for graphene charged with Q=+0.4 
e/cell. We also note that the excessive charging of single 
layer materials considered in this paper lead to instabil- 
ity. This is revealed from phonon dispersion calculations. 
For example, neutral graphene which has normally sta- 
ble planar structure and positive frequencies of acoustical 
branches in BZ, starts to attain imaginary frequencies of 
long wavelength acoustical modes at excessive charging. 
Weakening of graphene layer is expected to be reflected 
to its elastic properties, in particular to its stiffness. [39] 



to descend as a result of negative charging {Q < 0) and 
eventually they touch the Fermi level. Upon increasing 
Q these parabolic bands start to be occupied and hence 
part of Q is transferred to the quantum well in the vac- 
uum region. This way the rate of accumulation of ex- 
cess charge in the conduction band of single layer nanos- 
tructure recedes. Even if these parabolic bands appear 
as touching the Fermi level in the same band structure 
in (fc^:, fcj^) -plane they are physically separated from the 
states of single layer honeycomb structure under study. 
As it was mentioned before, this situation is an artifact 
of periodic boundary conditions and can be viewed as the 
vanishing of the work function. We note that in the case 
of negatively charged, finite size single-layer, the excess 
electrons are hindered from spilling out to the vacuum by 



IV. EXFOLIATION OF LAYERED BN AND 
MOS2 

We next investigate the exfoliation of single layer BN 
and M0S2 from their layered bulk crystal through charg- 
ing. We model 3-layer slab (sheet) of BN and M0S2 as 
part of their layered bulk crystal. We considered only 3- 
layer slabs in order to cut the computation time, since the 
model works also for thicker slabs consisting of 6-10 lay- 
ers graphene. '15] Energy minimizations of neutral sheets 
relative to stacking geometry are achieved. Stacking of 
3-layer BN and M0S2 slabs comply with the stacking of 
layers in 3D layered BN|3 and M0S2 crystals. [9j In these 
slabs, the layers are hold together mainly by attractive 
vdW interactions of a few hundreds meV and any repul- 
sive interaction overcoming it leads to exfoliation. When 
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Figure 4: (Color Online) (a) Variation of the ratio of lattice 
constants a of positively charged single layer graphene, BN, 
CH, CF and M0S2 to their neutral values ao with the average 
surface charge density, a. The unit cell and the lattice vectors 
are described by inset, (b) The charge density contour plots 
in a plane passing through a C-C bond, (c) Same as B-N 
bond. 



Table I: Dependence of the threshold charges on the vacuum 
spacing s (A) between 3-layer slabs. Threshold charge, Qe 
(e/cell) where exfoliation sets in and corresponding threshold 
average surface charge density = Q^jA (C/m^) are cal- 
culated for positive charged 3-layer Graphene, BN and M0S2 
sheets for s=50 A and s=20 A. The numbers of valence elec- 
trons per unit cell of the slab are also given in the second 
column. 



System 


#ofe 


Qe (e/ceU) 


G (C/m^) 


3-layer Graphene (s=50) 


24 


-H0.160 


+0.49 


3-layer Graphene (s=20) 


24 


+0.205 


+0.62 


3-layer BN (s=50) 


24 


+0.225 


+0.66 


3-layer BN (s=20) 


24 


+0.320 


+0.94 


3-layer M0S2 (s=50) 


54 


+0.322 


+0.57 


3-layer M0S2 (s=20) 


54 


+0.480 


+0.86 



electrons are injected to or removed from the slab, the 
Fermi level shifts up or down and cross the conduction or 
valance band of the insulator and attribute to it a metal- 
lic character. At the end, the excess charge by itself 
accumulates on the surfaces of the metallic slab inducing 
a repulsive Coulomb interaction between the outermost 
layers of the slab. Here we consider positive charging 
only, since in the case of negative charging the excess 
charges quickly spill into the vacuum before the exfolia- 
tion sets in. 

The amount of charge in the unit cell, which is neces- 
sary for the onset of exfoliation, is defined as the thresh- 
old charge Qe- Threshold charges are calculated for 3- 
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Figure 5: (Color Online) Variation of planarly averaged posi- 
tive excess charge X{z) along «-axis perpendicular to the BN- 
layers calculated for different s. As s increases more excess 
charge is transferred from center region to the surface planes. 
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Figure 6: (Color Online) Variation of cohesive energy and 
perpendicular force F±_ in 3-layer BN slab as a function of the 
distance L between the surfaces. Energies and forces are cal- 
culated for different levels of excess positive charge Q (e/unit 
cell). Zero of energy corresponds to the energy as i — > 00. 



layer slabs of graphene, BN and M0S2 for s=20 A and 
s=50 A. Results presented in Tab|l] indicate that the 
amount of threshold charge decreases with increasing s. 
This confirms our arguments in Sec. Ill that in positive 
charging large vacuum space, s, is favored. The mech- 
anism underlying this finding is summarized in Fig. [5] 
where we show the linear charge density, \{z) calculated 
for different s values of a 3-layer BN. For small s, the 
excess charge accumulates mainly at surfaces of the slab, 
also with some significant fraction inside the slab. How- 
ever, as s increases some part of Q is transferred from 
inside to the outer surface giving rise to the increase of 
the charge accumulation at the surface. At the end, for 
the same level of charging the induced Coulomb repulsion 
increases with increasing s. Accordingly, the same slab 
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requires relatively smaller amount of threshold charge Qe 
to be exfoliated, if s is taken large. 

In Fig. [6] we present the variation of the cohesive en- 
ergy of the 3-layer BN slab relative to three free BN layers 
for neutral Q — Q and positive charged Q > cases as 
a function of the distance L between the outermost BN 
atomic planes of 3-layer BN slab. The cohesive energy 
for a given L is obtained from the following expression: 
Ec = Et [3-Layer BN] - 3£'t [single layer BN]. The 
total energy of the single layer BN, Et [single layer] is 
calculated in a smaller supercell to keep the density of 
the background charge the same. The cohesive energy of 
the neutral slab in equilibrium is ^ 302 meV/cell. If the 
spacings of layers (i.e. L) starts to increase, an attrac- 
tive force Fj_ = —dEx/dL acts to restore the equilibrium 
spacing. F±{L) first increases with increasing L, passes 
through a maximum and then decays to zero. In Fig.|6] we 
also show how the minimum of cohesive energy decreases 
and moves to relatively large spacings with increasing Q. 
Concomitantly, the maximum of the attractive force for 
a given Q, F±^rnax decreases with increasing Q and even- 
tually becomes zero. This give rise to the exfoliation. We 
note that despite the limitations set by the neutralizing 
uniform charge on the total energy, the cohesive energies 
calculated for different charge levels reveal useful quali- 
tative information on the effects of charging. 

In Fig. [t] (a) we show isosurfaces of excess positive 
charge densities of 3-layer BN and M0S2 slabs. These 
slabs become metallic upon extracting electrons (i.e. 
upon positive charging) and excess charges reside at both 
surfaces of slabs. As shown in Fig.[7](b), the total energy 
raises with increasing charging or average charge density, 
a. In compliance with Fig. |6j the separation between 
surface layers, L, increases. The sharp drop of AE at 
or CTe indicate the onset of exfoliation due to the repul- 
sive Coulomb force pulling them to exfoliate. In Fig. [7] 
(c) L increases with increasing charging as discussed in 
Fig. [6] The increments of L exhibits a stepwise behavior 
for BN. This is also artifact of the method, where forces 
are calculated within preset threshold values. 

The variation of L of M0S2 slab with Q > display 
a different behavior due to charge transfer from Mo to 
S atoms. The exfoliation due to the static charging can 
be explained by a simple electrostatic model, where the 
outermost layers of slabs is modeled by uniformly charged 
planes, which yields repulsive interaction independent of 
their separation distance, i.e. F cx Q'^/{A-eQ), where cq 
is static dielectric constant. [T5^ Calculated forces differ 
from the classical force due to screening effects of excess 
charge residing inside the slabs. 

V. DISCUSSIONS AND CONCLUSIONS 

In this study, the threshold values of static charge, Qe, 
to be implemented in the slabs to achieve exfoliation are 
quite high. Such high static charging of layers can be 



achieved locally through the tip of Scanning Tunnelling 
Microscope or electrolytic gate. [38] The dissipation of lo- 
cally created excess charge in materials may involve a 
decay time td- Relatively longer to can induce a local in- 
stability and the desorption of atoms from nanoparticles. 
Experimentally ultra-fast graphene ablation was directly 
observed by means of electron crystallography. jl6| Car- 
riers excited by ultra-short laser pulse transfer energy to 
strongly coupled optical phonons. Graphite undergoes a 
contraction, which is subsequently followed by an expan- 
sion leading eventually to laser-driven ablation. [TQ Much 
recently, the understanding of photoexfoliation have been 
proposed, where exposure to femtosecond laser pulses 
has led to athermal exfoliation of intact graphenes. [TT] 
Based on time dependent DFT calculations (TD-DFT), 
it is proposed that the femtosecond laser pulse rapidly 
generates hot electron gas at ^ 20.000 K, while graphene 
layers are vibrationally cold. The hot electrons spill out, 
leaving behind a positively charged graphite slab. The 
charge deficiency accumulated at the top and bottom sur- 
faces lead to athermal excitation. [iTj The exfoliation in 
static charging described in Fig. [7] is in compliance with 
the understanding of photoexcitation revealed from pre- 
vious TD-DFT calculations, [ITj since the driving force 
which leads to the separation of graphenes from graphite 
is mainly related with electrostatic effects in both meth- 
ods. 

In summary, the present study investigated the effects 
of charging on the structural and electronic properties 
of single layer graphene, graphene derivatives, BN and 
M0S2, which have honeycomb structure. We concluded 
that while caution has to be exercised in the studies in- 
volving negative charging using large vacuum spacing, 
positive charging can be treated safely using large vac- 
uum spacing. 

We found that upon positive charging the band gaps of 
single layers of BN and M0S2 increase and the unit cells 
are enlarged. Consequently the phonons become softer. 
The charging of BN and M0S2 slabs were also studied. 
While these slabs are wide band semiconductors, they be- 
come metallic upon positive charging. Consequently, ex- 
cess charges are accumulated on the surfaces of slabs and 
induce repulsive force between outermost layers. With 
increasing positive charging the spacing between these 
layers increases, which eventually ends with exfoliation. 
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Figure 7: (Color Online) Exfoliation of outermost layers from layered BN and M0S2 slabs by positively charging of slabs, (a) 
Turquoise isosurfaces of excess positive charge density, (b) Change in total energy with excess surface charge density, (c) 
Variation of L of slabs with charging. 
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